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ABSTRACT 

We investigate the dynamics and radiation from a relativistic blast-wave which de- 
celerates as it sweeps up ambient matter. The bulk kinetic energy of the blast-wave shell 
is converted into internal energy by the process of accreting external matter. If it takes 
the form of non-thermal electrons and magnetic fields, then this internal energy will 
be emitted as synchrotron and synchrotron self-Compton radiation. We perform ana- 
lytic and numerical calculations for the deceleration and radiative processes and present 
time-resolved spectra throughout the evolution of the blast-wave. We also examine the 
dependence of the burst spectra and light curves on various parameters describing the 
magnetic field and non-thermal electron distributions. We find that for bursts such as 
GRB 910503, GRB 910601 and GRB 910814, the spectral shapes of the prompt gamma- 
ray emission at the peaks in vF v strongly constrain the magnetic fields in these bursts to 
be well below (<; 10~ 2 ) the equipartition values. These calculations are also considered 
in the context of the afterglow emission from the recently detected gamma-ray burst 
counterparts. 



1. Introduction 

The recent Beppo-SAX observations of fading X-ray afterglows coincident with the positions of 
gamma-ray bursts GRB 970228 (Costa et al. 1997) and GRB 970508 (Piro et al. 1998) have led to 
the first identifications of possible burst counterparts in the optical wave band (van Paradijs et al. 
1997; Djorgovski et al. 1997). In the optical spectrum associated with the latter burst, absorption 
lines have been detected with a redshift of z = 0.835 providing, for the first time, direct evidence 
for a gamma-ray burst at cosmological distances (Metzger et al. 1997). The gamma-ray fluence of 
GRB 970508 measured by BATSE is ~ 3 x 10" 6 erg cm" 2 (Kouveliotou et al. 1997). Therefore, 
if the absorption line measurements give a lower limit on the redshift of the burst, this implies an 
isotropic burst energy of > 10 51 erg s~ 1 . An impulsive event releasing this amount of energy in a 
compact region naturally leads to a fireball and thence to a relativistic blast-wave. 
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Blast-wave models for gamma-ray bursts have been examined previously in the literature in 
several contexts. The most extensive body of work on this topic has been produced by Meszaros, 
Rees and collaborators (e.g., Meszaros &; Rees 1992a; Rees & Meszaros 1992; Meszaros, Laguna 
& Rees 1993; Wijers, Rees & Meszaros 1997; Panaitescu k Meszaros 1998a, 1998b; et al.). Other 
recent papers include Sari & Piran (1995), Vietri (1997), Waxman (1997), and Katz & Piran (1997). 
The basic fireball/blast-wave model consists of some triggering event — either coalescing neutron 
stars or black holes (Meszaros & Rees 1992b) or the collapse of a massive star (Paczyhski 1998) 
or a failed type II supernova (Woosley 1993) — depositing a large amount of energy, Eq ~ 10 51 - 
10 55 ergs, in a small region with radius r$ ~ 10 6 -10 7 cm. Because of the unavoidable presence 
of baryonic material, it is expected that the initial fireball energy will be transformed into kinetic 
energy of these baryons rather than escape as radiation. This material expands until the internal 
motions of the baryons become sub-relativistic in the co-moving frame of the material, at which 
point it forms a cold shell with bulk Lorentz factor r ^ Eq/Mqc 2 , where M is the rest mass of 
the contaminating baryons. 

This shell continues to expand freely into the surrounding ambient medium until the integrated 
momentum impulse upon the shell by the swept-up matter is equal to the rest mass of the original 
material, ro47rr^p cxt /3 ~ Mq where p e xt is the mass density of the external medium. This defines 
the so-called deceleration radius (Rees & Meszaros 1992). Beyond this radius, the shell can no 
longer be regarded as freely expanding, and the bulk kinetic energy of the blast-wave begins to be 
reconverted into internal energy. If this internal energy is radiated promptly, then the deceleration 
of the blast-wave shell can be approximately described by T(r) oc r~ 3 , and the expansion is said 
to be in the radiative regime. On the other hand, if the internal energy is radiated on a time scale 
which is long compared to the expansion time scale, then the expansion is in the non-radiative 
regime, and T(r) cx r~ 3 / 2 . In either case, for large initial Lorentz factors, Tq ~ 10 2 -10 3 , relativistic 
effects (e.g., Rees 1966) compress the time scale for the radiation such that the bulk of the blast- 
wave energy is emitted in the first tens of seconds in the observer's frame following the initial 
detonation event, thus producing the observed gamma-ray burst. 

Several authors have pointed out that well after the prompt gamma-ray burst event the blast- 
wave shell will continue to decelerate and radiate (Vietri 1997; Waxman 1997). The recent detec- 
tions by the X-ray, optical and radio communities of the aforementioned fading X-ray, optical and 
radio counterparts within the error boxes of GRBs appear to support this picture. Furthermore, 
model estimates of the temporal decay of these transients yield time-dependencies which agree 
with those observed, F u ~ t^ 1 (Wijers et al. 1997); and for the one of the bursts for which optical 
data are available (GRB 970508), the optical spectral index is consistent with synchrotron emission 
from a power-law distribution of electrons, dN/dj oc with p ~ 2-2.3 (Djorgovski et al. 1997) 
indicating, for example, a shock-accelerated electron population. 

Despite the successes of the blast-wave model in accounting for the prompt burst properties 
and its prediction of fading afterglows, several important theoretical questions must be addressed 
in order to have a reasonably complete model: 
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• What is the nature of the coupling between the electrons, protons and magnetic field (Panai- 
tescu &; Meszaros 1998b)? To what extent can these components be in equipartition given 
that only the electrons can efficiently radiate away their energy? 

• How does the magnetic field change as the blast-wave decelerates? If it is initially formed 
through equipartition processes, but is not strongly coupled to the non-thermal electron 
energy density, how does it evolve? 

• What is the nature of the acceleration mechanism? Are the electrons energized by repeated 
diffusion across the shocks themselves or by gyroresonant scattering with disturbances in the 
post-shock turbulent MHD fluid? 

• What is the proper form for the injected electron energy distributions? Is it well described 
by a typical power-law spectrum? If so, what determines the characteristic energies of the 
particles? 



In order to begin to address these questions, it is worthwhile to go beyond the simple, though 
useful, back-of-the-envelope estimates which have generally prevailed in the literature thus far. In 
this paper, we present a detailed calculation of the blast-wave deceleration and the evolution of the 
magnetic fields and electron distributions under various assumptions relevant to the above issues. 
We compute model light curves and spectra and use the available afterglow data to help discriminate 
between the various options. The format of the paper is as follows: In § 2, we describe the basic 
dynamics of an impulsively driven blast-wave which decelerates by accretion of ambient material. 
In § 3, we discuss the physical processes responsible for producing the observed radiation including 
prescriptions for magnetic field generation, the formation of non-thermal particle distributions and 
the relevant radiation processes. The numerical procedures for computing the deceleration of the 
blast-wave and the integration of the blast-wave shell emission are described in § 4. Model spectra 
and light curves for gamma-ray bursts are presented and analyzed in § 5. Lastly, in § 6, we discuss 
these results, explore prospects for further research and present our conclusions. 



2. Deceleration of Impulsively Driven Blast-Waves 

We assume that a shell of material, with initial proper mass Mo, is injected with a bulk Lorentz 
factor To and that its subsequent deceleration is determined by energy and momentum conservation 
as it accretes mass in the form of matter and internal energy from the ambient medium and loses 
inertia via radiative processes such as synchrotron or inverse Compton emission. The internal 
energy of the material is assumed to be initially negligible compared to its rest mass energy. This 
corresponds to a time well after the random internal energy of the initial fireball has been converted 
into the bulk kinetic energy of the spherically expanding shell. We model the blast-wave as having 
a cross-sectional area which is a function of the distance traveled r. This functional dependence 
can be parametrized as a power-law: A(r) = Ao(r/ro) a = 47r/fer a where tq and Aq are the initial 
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radius and surface area of the blast-wave, respectively, and /& is a collimation factor to account 
for possible asymmetry of the explosion. In the calculations presented throughout this paper, we 
consider a spherical blast- wave with A{r) = Airr 2 . 

A straight-forward application of energy and momentum conservation gives the equation of 
motion of the blast- wave shell as it sweeps up ambient matter. We obtain the following expression 
(Blandford & McKee 1976): 

Ar r 2 — 1 

BL = -t t (i) 

dm M K ' 

where T = T(r) is the bulk Lorentz factor of the material, M = M(r) is the total mass including 
internal kinetic energy, dm(r) = A(r)p(r)dr is the rest mass swept-up in a distance dr, and the 
distance traveled r is measured in the stationary frame of the explosion. This frame will be referred 
to hereafter as the lab frame. 



2.1. Radiative and Non-radiative Limits 

As the blast-wave shell decelerates, its bulk kinetic energy is converted to internal kinetic 
energy. If this internal energy is not radiated away then it makes a significant contribution to the 
inertia of the blast-wave shell. In this case, the mass M, in equation ||, is given by 

M nr {r) = M + T Y{f)A{f)p{f)dr, (2) 

where the subscript nr identifies this as the non-radiative case. In the radiative case, the internal 
kinetic energy is immediately lost, and the inertia of the blast- wave shell consists only of the rest 
mass: 

M r (r) = M + f A{r)p{f)dr. (3) 

Note that since any radiation is assumed to be isotropic in the co-moving frame of the blast-wave 
shell, the bulk Lorentz factor is unaffected by this radiation and equation |l| still applies. 

In the radiative case, dM r = dm, the equation of motion (eq. [l]) is separable, and the resulting 
integrals are elementary. The solution is 

rr x = (M r (r)/M ) 2 (ro + l) + r -l 
[) (M r (r)/M ) 2 (r + l)-r + l u 

(Blandford & McKee 1976). In the non-radiative case, the deceleration is described by a pair of 
coupled ordinary differential equations: 

dF A(r)p(r)(T(r) 2 - 1) 

dr M nr (r) ' 1 ' 

dM nr 



dr 



T(r)A(r)p(r). (6) 
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These equations also have an analytic solution. The general expression for the mass is 



M nr (r) = M 2 + 2T M M sw {r) + M sw (rf 



1/2 



(7) 



where M sw (r) = J rQ drA(f)p(f) is the swept-up rest mass of the ambient material. The bulk Lorentz 
factor is given by 

M sw (r)+T M 
{ ] M nr {r) ■ W 

For a spherical blast-wave in a uniform external medium, A(r)p(r) = A7rr 2 pQ. For r ^> tq, the mass 
of the blast-wave shell is given by 

1/2 



M nr (r) 

and the bulk Lorentz factor is 



4tt 



r(r) 



1 



M nr (r) 



-yr 3 p + F M 



0) 



(10) 



From the above expression for the mass in the non-radiative case, we identify three regimes 
for the evolution of the blast-wave depending on which of the three terms dominates in equation 
For M > 2r M s ™, 



M nr {r) 
r(r) 



M , 
F . 



(11) 
(12) 



This regime corresponds to the initial period of free-expansion. For 2TqM sw > Mq > M sw /2Tq, 



M nr (r) w [2r M M^] 1/2 , 



r(r) 



M,. 



Mq y/ 2 1 

2T M SW ) T \2T M 



1/2" 



cx r 



-3/2 



(13) 
(14) 

(15) 



Here we recover the familiar power-law r _3//2 dependence of the bulk Lorentz factor which is often 
quoted for a non-radiative relativistic blast-wave decelerating in a uniform medium. Finally, when 
M sw > 2FqMo the bulk motion of the blast-wave is non-relativistic, and we have 



M nr (r) w M sw (r), 



r(r) 



i + r c 



■Mo 

M sw (r) 



(16) 
(17) 



The velocity of the blast-wave inferred in this limit, using the approximation T ~ 1 + /3 2 /2, is 
/3 = \/YqMq/M sw cx r -3 / 2 , which agrees with the adiabatic solution for non-relativistic blast-waves 
(Lozinskaya 1992). 
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Energy conservation gives the rate of accreted kinetic energy as a function of time in the lab 
frame: 

— = c 3 A(r)p(r)P(F 2 -T). (18) 

This expression follows directly from the equation of motion (eq. |j) and therefore applies regardless 
of whether the blast-wave is radiative or non-radiative. Since the above quantity is a Lorentz 
invariant, it also gives the injected power in the co-moving frame of the blast-wave (Blandford & 
McKee 1976). 



2.2. Partially Radiative Regimes 

In general, the evolution of the blast-wave shell will be intermediate between the radiative and 
non-radiative limits. If we denote the radiated power in the co-moving frame as dE Ta ^/dt, then the 
relativistic mass of the material is 

M(r) = M + f r drT(r)p(r)A(f) - f dt\^±, (19) 

Jr JO C dt 

where the latter integral is performed over co-moving time and represents the relativistic mass lost 
to radiation. Converting to the lab frame coordinate using dr = c/3Tdt, and differentiating wrt r, 
we obtain a more general form for equation ||: 

This expression must be integrated numerically in conjunction with the equation of motion (eq. [I]). 
The radiated power dE vax i/dt is determined by the co-moving particle distributions and the relevant 
energy loss mechanisms. We address these issues in the following sections. 



3. Physical Processes 

3.1. Particle Injection Distributions 

We assume that some fraction, £ e , of the accreted kinetic energy is injected directly into 
non-thermal electrons. Energization of these electrons may occur via acceleration by shocks (see 
e.g., Ellison, Baring & Jones 1995) or by gyroresonant scattering with plasma- wave turbulence 
induced by the swept-up material (Dermer, Miller & Li 1996 ). However, since we do not model 
these processes in detail, we adopt a power-law injection spectrum for the electron momenta, 
dN/dp ~ p~ s . In terms of the initial electron Lorentz factors, 7j, this injection spectrum per unit 
co-moving time is 

"' V d %( 7 2 _ l)( S +D/2. (21) 



djdti dti 
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Here the time-dependent normalization of the injection spectrum is given by the prescribed fraction 
of the power available as kinetic energy: 

dN _ Z e dE/dt 



dt 



/7't,max 
(7 2 -7*)(7*-l) (s+1)/2 ^ 
,'i,min 



The minimum injected electron Lorentz factor in the co-moving frame should be at least as large as 
the bulk Lorentz factor of the blast -wave, 7^, mm — P* However, if the electrons are in ecjuipartition 
with the protons, then the lower cut-off for the electron distribution will be approximately the 
typical energy of the incident protons in the co-moving frame, 7i jm i n = (m p /m e )T (Katz & Piran 
1997; Panaitescu Sz Meszaros 1998a). We parametrize the uncertainty between these two extremes 
with a parameter r\: 

Tvnin = T/^rfa). (23) 

The upper cut-off of the electron distribution is given by balancing the radiative losses with 
the power gained due to the acceleration. By balancing acceleration time scales for both shock and 
gyroresonant processes with synchrotron loss time scales, de Jager et al. (1996; see also de Jager 
& Harding 1992) argue that for a given magnetic field the maximum electron Lorentz factor is: 



7i 



~ 4 x W 7 (B/1 Gy 1/2 . (24) 



3.2. Magnetic Field Prescriptions 

At a minimum, the post-shock magnetic field will be that of the ambient external field, am- 
plified by the compression ratio (~ 4 for a strong shock) and the Lorentz transformation to the 
co-moving frame: B ~ 4n? pre ~ 10~ 5 r Gauss. Here we have assumed an ambient magnetic field 
intensity typical for the ISM of a few fiG. Even for bulk Lorentz factors of T ~ 10 3 , this field 
strength is too low for the electrons to be radiatively efficient. Another source of the magnetic field 
energy is the remnant field which was created during the initial fireball event and which gets carried 
along with the expanding blast-wave material. However, if flux freezing holds, i.e., r 2 B ~ constant, 
the field energy will decrease in proportion with the internal energy of the contaminating baryons, 
and by assumption we only consider the shell after it has expanded to the point at which the shell 
material is cold in the co-moving frame, hence the remnant field energy density should also be 
negligible. 

If gyroresonant processes mediate the approach to equipartition between the protons and 
electrons, then it follows that the magnetic field will also be in equipartition with the particles. 
In true equipartition, the bulk of the internal energy will be deposited in components which are 
themselves not efficient radiators, namely the protons and the magnetic field. However, if strong 
coupling exists between the electrons, protons and the magnetic field such that these components 
remain in equipartition even as the electrons efficiently radiate away their kinetic energy, then we 
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would expect the dynamics of the blast-wave to be in the radiative regime (Panaitescu & Meszaros 
1998b). In practice, this would correspond to £ e ~ 1 since the electrons would then be able to tap 
all of the available internal energy of the blast- wave shell. If the coupling between the magnetic field 
and particles is weak, then we would expect the internal energy in the protons and magnetic field 
to remain in the blast-wave and the blast-wave dynamics should be intermediate between radiative 
and non-radiative. In either case, we would expect that the generated magnetic field energy density 
will be proportional to the energy density of the swept-up material. This implies a magnetic field 
of the form: 

B= (327r£ B Pext) 1/2 Tc, (25) 

where we have assumed a compression ratio of 4. The parameter contains all the uncertainties 
regarding particle-field coupling and the degree to which equipartition is achieved and maintained. 
In principle, £b will change with time depending on the complicated plasma physics processes taking 
place within the blast-wave. However, we will treat it as a constant and address the implications 
of this simplification in our discussion below. 



3.3. Radiation Processes 



The radiation processes which we consider are synchrotron and synchrotron self-Compton. An 
expression for optically thin emission spectra for these processes are given by Dermer, Sturner & 
Schlickeiser (1997) where they use a 5- function approximation for the synchrotron spectrum of a 
single electron which is centered at the cyclotron frequency and which is normalized by the total 
synchrotron power. When convolved with a smooth distribution of electron energies, this approxi- 
mation provides a reasonably accurate representation of the optically thin synchrotron emissivity: 

dE 



isyn(e) 



dtdVdedn 
co T u B /_e_\ 
6vre B \e B ) 



V2 dN 



dVdj 



e 



1/2 



(26) 



Here e = hu/m e c 2 , €b = B/B cvlt = 5/4.414 x 10 13 G, <jt is the Thomson cross-section, and 
dN/dVd^f is the electron distribution function. Unfortunately, this approximation breaks down at 
the end-points of the electron distribution. In order to obtain a more accurate photon spectrum, 
we apply a more sophisticated calculation (Crusius & Schlickeiser 1986; 1988) appropriate for an 
isotropic distribution of electrons. However, we note that since the sharp electron cut-offs are 
themselves approximations, the spectra and light curves due to electrons at these cut-offs should 
be regarded with caution. 

Using a similar (5-function approximation, Dermer et al. (1997) also derive an expression for 
the first-order SSC emissivity in the Thomson regime: 

rmin(e,l/e) 



Jssc(e) 



2 1/2 

ca^UBr^e 1 



9vre' 



3/2 



Je B 



die' 



1 dN 




dN 




dVdj 




dVd-f 





(27) 
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The factor r& comes from the characteristic escape time (~ rjj/c) for synchrotron photons to leave 
a homogeneous spherical blob. In the case of a blast-wave shell we take to be the thickness of 
the shell in the co-moving frame. 



4. Numerical Procedures 

4.1. Evolution of the Non-Thermal Particle Distributions 

The non-thermal particle distribution is governed by the equation of continuity in energy space: 

d fdN\ d f.dN\ dN 

The losses consist of radiative losses due to synchrotron and synchrotron self-Compton emission 
and adiabatic losses: 

7 = 7rad + 7adi (29) 

= -^2^VccT T ('UB+'u syn (7))- |— . (30) 

The first term on the rhs corresponds to the radiative losses where u syn is the energy density of 
the synchrotron photons which can be up-scattered by an electron with Lorentz factor 7 in the 
Thomson limit: 

° l h , , dn 
dedt 

Here dn syri /dedt is the synchrotron photon number spectral density. The second term on the rhs of 
equation |3(] is the expansion losses where V is the volume of the emitting material. However, since 
the relevant radiative loss time scales are much shorter than the expansion time scales, we neglect 
the adiabatic loss term. We also find that bremsstrahlung losses can be neglected compared to the 
synchrotron and SSC losses. 



^syn(7) = / m e c 2 dee s J n . (31) 
Jo 



The evolution of the particle distribution is computed by integrating equation 28 using an 
implicit finite differencing scheme (Press et al. 1992): 

N i+ ij - Njj + j i+lj+1 N i+lj+1 - 'y i+lj N i+lj = ^ ^ 
U+i - U 7j+i - 7j J 

Here Nij = dN/dj(ti,jj), jij = 7(^,7^), and the source term, Qij, is the injection function given 



by equation 21. This implicit finite differencing scheme yields a tridiagonal system of equations 
which is readily inverted to give the electron distribution as a function of time. We note that the 
losses associated with the SSC term depend on the electron distribution function through n syn and 
must be determined self-consistently. This is achieved by iteration until the electron distribution 
has converged. 
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4.2. Self- Consistent Dynamics 



If the kinetic energy deposited in the electrons is immediately radiated, then equation 2C is 
equivalent to 

AM 

rtr)A(r)[(r(r)-l)(l-£e) + l]. (33) 



dM 
dr- 



ill general, however, dE Ta ^/dt depends on all previous history of the accretion process through 
dN/d'y and the integration of the equation of continuity (eq. p^), hence dE vaj ^/dt contains an 



implicit integral over r. We integrate equation |20|by taking time steps equal to the shortest energy 
loss time scale At ~ (7/7) m in> and for each of these steps, estimate dE raf ±/dt using 



rf-Erad \ __ m e c 2 
dt ) t .~ (At)i 



c?7 7 



dN 
dj 



(-) 



(34) 



where t{ = U-\ + At is the time of the ith step. The bulk Lorentz factor and total mass are 
then obtained by forward-stepping: Tj = Yi-\ + Ar(dT/dr) and Mi = Mj_i + Ar(dM/dr), where 
Ar = c/3TAt. The estimate of the radiated power (eq. U) from this simple forward-stepped 
integration is then used as input for a variable step-size Runge-Kutta integrator and the process is 
repeated until convergence. Only for very low magnetic fields, <J of the equipartition value, 
are several additional iterations required in order to obtain convergence. 



4.3. Computing the Observed Spectra from the Co-Moving Emissivities 

Using the electron distribution function from the integration of the equation of continuity 
and the expressions for the synchrotron and SSC emissivities, we compute the radiation in the 
co-moving frame of the blast-wave shell. Because of light travel time effects, emission from the 
blast-wave shell at different lab frame times contributes to the observed flux at a given observer 
time. Hence, in the observer frame, the flux at any given time will consist of an integration of the 
emission over the course of the blast-wave expansion. If a shell element moves radially outward 
from the point of the explosion with velocity c/3(r) and at an angle 9 = cos -1 \i with respect to 
the observer line-of-sight, then the time delay between radiation from this shell element at radius 
r and the initial explosion event is given by 



^obs 

C 



df 

W) 



r cos ( 



(35) 



The observed photon flux is 
/ dN \ 



\dAdedtdQ J 



obs 



dr 



dA ^ ( lA^uin ) 6 ( r - f (*obs,^)) 
shell \dAdedtdtt J co 



de c 
de 



obs 



dt co 

^obs 



On 



dA C0 
~dA 



(36) 
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where di is the luminosity distance to the blast-wave, the cL4-integral is performed over the shell 
surface at radius r, the element of area is dA = 2irr 2 dp, and the factor of \i accounts for the 
projection of the surface area of the shell. Note that for convenience we have taken the redshift to 
be z = 0. It is straight-forward to include the appropriate number of (1 + z) factors in the final 
expressions. The r in the 5- function satisfies the time delay equation (^5[) for given values of i bs 
and fj,. The transformation from the co-moving frame to the observer frame introduces a number of 
Jacobian factors, all of which are powers of the Doppler factor T> = [T(l — \d£l co /dQ\ = T> 2 , 

\de bs/de co \ = V, \dt Q ^ s /dt co \ = V^ 1 , and \dA co /dA\ = 1. The area is unchanged since the velocity 
of the shell is perpendicular to its surface. Explicitly including all these factors and using the fact 
that the emission is isotropic in the co-moving frame, we obtain 



/ dN \ If 1 ! 



\dAdedtdQ J obs Aird 2 J-i 8tt 



dN fj,d(j, 
de co dt co [T(l - (3n)f 



(37) 

=K*obs,M) 



where dN / de co dt co is calculated from the expressions for the synchrotron (eq. [2fJ) and synchrotron 



self-Compton (eq. |27| ) emissivities and the integration of the electron equation of continuity (eqs. 28 



At a time delay t Q b s the apparent size of the afterglow image on the sky will be given by the 
maximum value of the quantity p = ryl — p? . Differentiating wrt p and using the expression for 
constant time delay (eq. |35|) this maximum occurs at = (j,, so that p max = f/Y(f). The quantity 
f is found by solving the equation 

r df 

= ct obs- (38) 

Jo p(r) 

For constant (3, the solution is f = cY 2 j3t Q \ iS yielding the standard result p max = cY[3t \> s . For a 
blast-wave shell in a uniform external medium whose Lorentz factor can be described by V = Tq 
for r <rd = (3Mo/87rro/9o) ) an d T = Tq(t /V^) - ^ for r > r^, where C ~ 3/2 in the non-radiative 
case and ( ~ 3 in the radiative case, as long as T > 1 and i Q bs > (2C/(2C + l))*obs,o where 
^obs,o = r d/2rgC, the apparent size is 



fmax — p 
1 







1 '( 2C + l)i2^L_2C 



2(C + 1)\ *obs,0 



(2C+1)/(C+1) 

(39) 



In the radiative and non-radiative cases, (2( + 1)/(C + 1) = 4/7 and 5/8 respectively (cf. Waxman, 
Kulkarni & Frail 1997). For purposes of illustration, we take an intermediate value, (2£+l) / (C+l) = 
3/5, and obtain 

p m „_l.4xiu ^ 103 j ^10=1 ergs J Ucm-sJ { lsec ) cm ' W 

For the radio afterglow associated with GRB 970508, Frail et al. (1997) find that the decline of the 
8.64 GHz flux variability due to interstellar scintillation implies a source size of p ma , x ~ 10 17 cm at 



tobs ~ 1 month. Using equation 40 above, this observation is consistent with a spherical blast-wave 
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with initial energy Eq ~ 10 51 ergs expanding in an external medium with density n ex t ~ 1 cm -3 . 
However, the very weak dependence of p max on Tq, Eq, and n ex t means that any of these parameters 
could vary by as much as two orders of magnitude and not substantially change this result. The 
significance of equations [39] & ^ is the dependence of the apparent size with time. For a sufficiently 
bright radio afterglow, the blast-wave hypothesis could be tested by measuring the time dependence 
of the amplitude of the radio scintillation to determine if it conforms with the above prediction for 
the apparent size as a function of time (Goodman 1997; Frail et al. 1997; Waxman et al. 1997). 

5. Model Light Curves and Spectra 

We test our implementation of the above calculations by comparing our computer program's 
output with the analytic results presented by Dermer & Chiang (1998) for the case of a radiative 
blast-wave and a constant magnetic field. For these calculations, we use Tq = 300, B = 1 G, 
an injection index of s = 2, and the blast-wave is assumed to be radiative. Figure |l]a shows the 
spectra, e 2 dN/dedt, as seen by an observer at times t b s = 1, 10, 10 2 , . . . seconds for the numerical 
(solid curves) and analytic (dashed) calculations. The vertical lines indicate radio (2.4 GHz, dot- 
ted), optical (1 eV, solid), X-ray (1 keV, dashed) and gamma-ray (1 MeV, dot-dashed) energies. 
In figure |l|b, the light curves corresponding to these energies are plotted, with the results from the 
analytic calculation over-plotted as symbols. The agreement is reasonably good, and the discrep- 
ancies are due to the fact that the analytic calculations use T(r) = To(r /r^)" 3 for the bulk Lorentz 
factor rather than equation |], somewhat different expressions for the electron energy losses due 
to synchrotron radiation, the electron distribution function and the ^-function approximation for 
the synchrotron emissivity. In addition, the use of the finite-difference method for calculating the 
evolution of the electron energies results in less sharp features (e.g., rounded high energy cut-offs) 
in the distributions. 

In the upper panel of figure ^, we plot the surface brightness in the radio waveband as a 
function of perpendicular distance, p, for this calculation. As before, a curve is shown for time 
delays i b s = 1, 10, 10 2 , . . . seconds (plotted top to bottom). Each curve is computed as a function 
of observer angle 9 = cos -1 fj,, thus the surface brightness is double-valued and turns around at 
Pmax- The flat portion within p max results from emission from the leading surface of the blast-wave 
shell, i.e., for 9 < T(f). In the lower panel, we plot p ma . x for each time delay, i bs> along with the 
curve given by equation |3^ using Eq = 10 50 ergs, Tq = 300, n ex t = 1 cm" 3 , and £ = 3. 

An example of a more realistic calculation is shown in figure ||. In figure ^a, we plot as the 
solid curve the bulk Lorentz factor of the blast-wave shell as a function of distance from the central 
explosion. The initial Lorentz factor is Tq = 300, the fireball energy is Eq = 10 52 ergs, and the 
ambient density is n ex t = 1 cm -3 . We use = 1 and r/ = 1 for the minimum electron energy 
(eq. ^3|), and the index of the electron momentum distribution is s = 3. The radius at which 
the free-expansion phase ends and the deceleration phase begins is consistent with the analytic 
estimate of = (3Eo/87nripC i n ex tT'Q) = 2 x 10 16 cm. The fraction of kinetic energy made available 
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for radiation by electrons is £ e = 0.5. Since this is less than unity, the evolution of the bulk 
Lorentz factor lies between the two limiting cases with T(r) oc r~ 19 for r > r^. Figure ||b shows 
the photon spectra, including the synchrotron self-Compton component, at observed time delays 
tobs = 1) 10, 10 2 . . . seconds after the initial fireball event. The blast-wave radii corresponding to 
these times are indicated as dotted vertical lines in figure |3|a. 

At early times, the shape of the synchrotron part of these spectra can be inferred directly from 
the corresponding co- moving electron distributions at the leading surface of the blast-wave (i.e., at 
H = 1) since the bulk of the observed radiation originates from within a small angle ~ 1/T around 
this region. These distributions are shown in figure |||c with the earlier distributions corresponding 
to the lower curves. The shapes of these distributions are determined by the momentum dependence 
of the particle energy loss rates and the shape of the injection spectra. The steep sections at high 
energies, above the break but below the high energy cut-off, are the power-law-injected electrons 
which have evolved to a cooled distribution with dN/d'j ~ (Dermer & Chiang 1998). 

The break energy (e.g., 7t>reak — 6 x 10 5 for the t Q ^ s = 1 s curve) is the minimum injection 
energy corresponding to the leading part of the blast- wave, 7b re ak ~ 7i,min = r/(m p /m e )r, and the 
electrons with the smallest energies (7 m i n — 10 2 ) are those which were injected at early times with 
7i = (m p /m e )To but which have cooled due to radiative losses. The flat sections at energies less 
than 7 m i n result from synchrotron cooling of continuously injected electron distributions which have 
sharp low energy cut-offs. The spectral shape, dN/d'j ~ 7~ 2 , is due to the energy dependence of the 
energy loss rate, 7 oc — 7 2 . This shape is independent of the injection power-law index and hence 
will be a generic feature of a synchrotron-cooled electron spectrum resulting from injection spectra 
with sharp lower cut-off injection energies. An analytic expression for this part of the spectrum can 
be obtained by treating the injection function as a 5-function at the minimum injection energy: 

— Sl-yi-ri—^Ti) / djai , (41) 



d^idti dti \ m . 

where Tj = T(rj), n is the radius at which this distribution is injected and the integral on the rhs 
sets the normalization. We find, in terms of observing time t \, s , 



2 dN _ ( 3irm e ^ e \ J (i bsAobs,o) *obs < t ohsfi 

7 lh ~ \4a T m p Z BV )r~o X \ [(2C + l)(W*ob Sj0 ) " 2fl {C+2)/(2C+1) t ohs > t obSj0 . ' " ' 

This expression is valid for s > 2. The electron spectra implied by this expression, including the 
cooled high energy power-law distribution above the breaks at 7i,min and the estimated roll-overs 
at low energies, are plotted as dashed curves in figure [3|c. 

In figure |b, the energy of the peak of the synchrotron spectrum corresponds to the break at 
high electron energies. These quantities are related by 

B 2 

Cpeak ~ "5 7break r ( 43 ) 



crit 



327rf Rmr,c 2 n cxt r 2 \2 

' ' v 2 (^) r 3 (44) 



B„\ t V m. 
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« 3xlO- 8 ^ t Vr 4 (45) 

(Katz & Piran 1997). Using this relation in conjunction with the expressions for the synchrotron 
spectrum (eq. |26| ) and the electron distributions (eq. ^2|), we find the temporal dependence of the 
peak of the synchrotron component: 

2 dN \ 2 2 2 J (^obs/*obs,o) 2 ^obs < ^obs,0 

6 a^J pcak " 27FmpC X \ [(2C + l)(WU s ,o) - 2C] (2 " 2 « /(2?+1) t obs > Us.o- 

(46) 

Here we see that the maximum of the light curve of the peak energy occurs at i G b s o = r d/2cY'^. We 
also note that the above expression is independent of the magnetic field equipartition parameter. 

With an injection index of s yielding a power-law index for the cooled distribution of p = s + 1, 
the energy index of the corresponding synchrotron emission is a = s/2 (L v oc v~ a ). Depending on 
the strength of the magnetic field, the energy index of the synchrotron spectrum just below 6 P eak 
will be either a = 1/2 or —1/3. If the field strength is well below equipartition, i.e., £b 1, then 
the cooling time will be very long, the electrons at the low energy cut-off will have little time to cool, 
the flat portion of the spectrum will be very narrow or non-existent, and the low energy spectrum 
will be that of a cut-off electron distribution: j S yn cx e 1 / 3 (Katz 1994; Rybicki & Lightman 1979). 
This latter behavior is seen in the spectra of the bright bursts GRB 910503 and GRB 910814, while 
a softer transition region with j syn cx e -1 / 2 between the harder a = —1/3 section at lower energies 
and the peak in vF v may be present in the spectrum of GRB 910601 (Schaefer et al. 1998). This 
may indicate that the synchrotron cooling of this latter burst is relatively more efficient than that 
of the other two bursts. 

More quantitatively, the location of the spectral peak and the shape of the spectrum above and 
below e pea k place strong limits on the equipartition parameter the initial bulk Lorentz factor 
Tq and the injection spectral index s. For t Q ^ s > £ o bs,0) the lower cut-off electron Lorentz factor is 

-l 



7min = 70 



i _i_ a o*o h ( r X 



(47) 



where ao = (16/3)caT(m p /m e )n ex t^B^o and 70 = 7](m p /m e )TQ. In deriving this expression, we 
assume that the particle injection is only significant once the blast-wave reaches the deceleration 
phase r > r^. For a burst spectrum in which electron cooling is important at all electron energies, 
which occurs if t co ~ 2Fo£ bs < (7/7)min> the minimum electron Lorentz factor obeys 

TYL 

7min < 7i,min = T) — T (48) 

m e 

= 7o (£) . (49) 



Writingr/r^ = l+(r— r^j/rd and expanding equations |47| fc [49| to first order, we obtain ao ^ C/7o*o- 
In terms of the relevant model parameters, 

V^rt /3 ^0.09E; 2 1/3 nf /3 C, (50) 
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where £"52 = E/10 52 erg and n% = n ext /lcm 3 . Using this expression in conjunction with equa- 
tion we find 



Thus, for burst spectra peaking at ~ 1 MeV and for which electron cooling is important, we obtain 
an upper limit for the initial bulk Lorentz factor of ~ 300. Remarkably, this upper limit depends 
only very weakly on the total burst energy and the density of the ambient medium. For bulk 
Lorentz factors less than the above limit, the equipartition fields will be sufficiently strong so that 
the energy spectral index spectral index just below e pea k will be a = 1/2 and there will be a break 
at even lower energies below which a = —1/3. 

Bulk Lorentz factors larger than the above limit are possible, but in these cases, the electrons 
will not be fully cooled. If the injection index is s > 3, then just above e pea k the energy spectral 
index will be that of the uncooled power-law electron distribution, a = (s — l)/2. Furthermore, if 
the upper injection cut-off, 7j,max> is sufficiently high, an additional spectral break of Aa = 1/2 will 
exist at higher energies marking the transition between the cooled and uncooled electrons (Dermer 
& Chiang 1998). Such a break would manifest itself as curvature in the high energy tail above 
epeak- If the electron injection index is s < 3, then the peak in vL v may actually correspond to the 
transition between cooled and uncooled power-law electron distributions. More likely, however, the 
peak in vL v will correspond to the synchrotron emission from the electrons at the upper cut-off 
(unless this upper cut-off energy is very high). In this case, there will be a spectral break between 
a = —1/3 and a = (s — l)/2 below the peak in uL u and an exponential cut-off above the peak 
(Rybicki & Lightman 1979). None of the bursts for which good spectra are available show an 
exponential cut-off above the vL v peak or a spectral index change of Aa = 1/2 at the peak itself 
(Tavani 1996; Schaefer et al. 1998). In figure |||, we show burst spectra illustrating the various 
spectral shapes we have described. 

The evolution of the electron distributions will also be reflected in the afterglow emission. 
In figure [3|d, we show the light curves at radio, optical, X-ray, gamma-ray, and TeV (dot-dot-dot- 
dashed curve) energies. As with our earlier calculation, these energies are indicated by vertical lines 
in the figure ||b. As the low energy electron cut-off evolves downward due to cooling, the hard, low 
energy part of the synchrotron spectrum will also move downward in energy. The evolution of the 
bulk Lorentz factor as the blast-wave decelerates contributes to this effect as well. At lower observer 
energies, the light curve will peak later and the spectral index will get softer with time. However, 
the precise shape of the light curve in any given band is difficult to predict due to the uncertainties 
in the nature of the electron cut-off and the evolution of the magnetic field. The general trend of 
hard-to-soft evolution should nevertheless hold, even if only qualitatively. Once the blast-wave has 
evolved so that the emission of the cooled power-law electron distribution lies within the optical 
band-pass, a measurement of the electron distribution power-law index can be inferred from optical 
afterglow data. For the afterglow of GRB 970508, the optical energy spectral index is a ~ 0.7 ±0.3 



r <s 



2.6 xio^V^V^V 9720 ^ 
s^xio-^^V^C^V^Sc, 




(51) 
(52) 



-16- 



(Djorgovski et al. 1997), which is consistent with a cooled power-law distribution with s = 1.4. 
This value is quite different from the injection indices implied by the prompt gamma-ray emission 
of GRBs 910503, 910601 & 910814 (Schaefer et al. 1998). The energy spectral indices for the high 
energy tail of the prompt burst emission range from a = 1.3-2.5, corresponding to injection indices 
of s = 2.6-5 if these electrons are cooled and even larger values of s = 3.6-6 if they are not cooled 
(Tavani 1996). 

In addition to the evolution of the electron distribution injection index, it is expected that the 
magnetic field will not stay at a constant fraction of its equipartition value. If there is some sort 
of time scale, T eq , for particle- field equipartition to obtain, then one might expect that £b ~ i~ eq /to 
where to = r d/cTo is the dynamical time scale of the blast- wave expansion in the co- moving 

-4 



frame. From equation 45, bursts which peak at 1 MeV have weak magnetic fields with £g ~ 10" 



for Tq = 300. This implies a rather long equipartition time scale of r eq ~ 2 x 10 7 s. Even 
for a non-radiative burst, this time scale is reached only after the burst has almost completely 
decelerated. Most likely, particle-field equilibration is a highly non-linear process and thus cannot 
be parametrized as a simple linear function of time. The lack of a proper theory describing such 
processes poses a major uncertainty in using the characteristics of the prompt gamma-ray burst to 
set model parameters for extrapolation of the behavior of the blast-wave to describe the emission 
at much later times. 

Nonetheless, for purposes of illustration, we present a calculation which largely matches the 
observed properties of typical prompt GRB emission. In as much as the spectral slopes for GRBs 
below the vF v peak are consistent with a. = —1/3 (Schaefer et al. 1998), the limiting values of 
ro = 260 and £b = 5.5 x 10~ 5 apply. These values also insure that for the prompt burst spectrum 
we have e pea k ~ 1 MeV. We have also set Eq = 10 52 ergs, and used a compromise value for the 
injection index of s = 3. In figure [5], we show the results of a calculation for those parameters. The 
blast- wave is almost non-radiative, with ( ~ 1.6. Despite the aforementioned uncertainties in the 
model, this calculation displays many of the properties which have been seen in actual bursts and 
their putative afterglows. In comparison with the calculation shown in fig. ||L it demonstrates the 
variety of afterglow behaviors which can be obtained. 

In figure |b, the gamma-ray emission in the prompt burst extends up to 1 GeV, well into the 
EGRET energy range (cf. Dingus 1995), the initial rises in the X-ray and optical light curves follow 
that of the prompt burst itself by ~ 10 3 s and ~ 5 x 10 4 s, respectively. However, because of the 
steeper injection index implied by the prompt burst, the flux in the X-ray band decays as ~ t , 
significantly faster than that measured for GRB 970508 (Piro et al. 1998). The optical light curve 
decays somewhat slower as ~ t 3 since the optical emission corresponds to the flatter part of the 
spectrum near the peak in vh v . The prompt burst spectrum evolves from hard to soft, and its 
time scale is ~ 10 seconds. The radio light curve is rising as ~ i 3 and the radio spectral index is 
a = —1/3. Furthermore, this model makes the prediction of TeV emission due to SSC processes 
which should peak at a time comparable to the peak in the X-rays. 
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6. Conclusions 

In this paper, we have attempted a more realistic calculation of the dynamics and synchrotron 
and synchrotron self-Compton emission for the blast-wave model of gamma-ray bursts. By matching 



the detailed characteristics of burst spectra, we have found relations (eqs. |5l| & 52) which place 
constraints on magnetic field strengths and bulk Lorentz factors. If these relations are to be 
believed, then burst data can have a significant impact on models of magnetic field generation in 
turbulent plasmas. 

The apparent deficiencies of this calculation point towards areas of further research. In partic- 
ular, the detailed temporal structure of individual bursts is not explicitly dealt with in this model. 
In the context of external shocks, it may be due to inhomogeneities in the external medium, or 
fluctuations in the electron injection and/or magnetic field equipartition parameters (Panaitescu 
& Meszaros 1998a). For s > 3, the generalized expression for the luminosity for energies e > e pea k 
(eq. |4|) is 

= 27rm p c 2 e e r 2 n ext (r)r 2 (r)(e/e P eak) A (53) 

where A = (2 — s)/2 applies for relatively strong magnetic fields £b <; 1CP 4 when the electrons 
just above the break are efficiently cooled and A = (3 — s)/2 applies for relatively weak fields and 
uncooled electrons. From this expression we see that any burst light curve substructure must be 
due to variations in £ e and n ext , and indirectly, due to variations in £b through e pea k (eq. |45| ). 

Our treatment of the dynamics also ignores the structure of the shock region itself, and our 
approach essentially only considers the emission from the forward shock and neglects the reverse 
shock. Panaitescu & Meszaros (1998a) have performed calculations similar to our own, but from a 
hydrodynamical perspective, and found that the reverse shock only makes a significant contribution 
to the emission at optical and UV energies. Therefore, neglecting the reverse shock should not affect 
our results for the gamma-ray emission, but it could have a significant impact on the optical and 
radio afterglow emission. We also neglect the thickness, Ar, of the shock shell in integrating 
the emission for a given observer time t b s - This should not be important at early times when 
Ar = to/Tq (in the lab frame), but it could affect the afterglow emission at late times. 

It is unlikely that the blast-wave itself is spherical. If the initial fireball is created by the 
coalescence of two compact objects, then the orbital plane defines a natural axis of symmetry along 
which the blast-wave will propagate (Meszaros & Rees 1992b). This sort of asymmetry could be 
accounted for in our model by a non-unity collimation factor, fa. Furthermore, if the observer 
line-of-sight does not lie within the opening angle of the blast-wave cone, then other effects due to 
relativistic beaming and the gradual deceleration of the shock front would be introduced. In this 
respect, highly anisotropic blast-waves would share properties with relativistic jets in blazars. 

This analogy can be take even further by noting the similarity of the burst spectra we derive 
compared to that of gamma-ray blazars. Like our model spectra, the spectral energy distributions 
(SEDs) of these objects tend to have two peaks, one in the UV-X-ray range and one at gamma-ray 
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energies. If the lower peak in blazar SEDs is due to synchrotron emission and corresponds to the 
~ 1 MeV peak in gamma-ray bursts, we can apply a similar analysis as we have discussed above 
to derive bulk Lorentz factors and equipartition parameters for blazars. In particular, the recent 
ASCA observations of Mrk 421 (Takahashi et al. 1996) provide sufficient information to get actual 
values rather than simply upper or lower limits. Using the light curves of Mrk 421 measured in 
different X-ray energy bands, Takahashi et al. (1996) performed a cross-correlation analysis and 
found that the longer relative time lags of the lower energy data versus the higher energy data are 
consistent with synchrotron cooling of the underlying electron distribution. Several authors have 
noted this effect and have calculated this temporal dependence for the cases of bursts and blazars 
(e.g., Tashiro et al. 1995; Tavani 1996; Dermer 1998). Takahashi et al. use the TeV variability time 
scale (Kerrick et al. 1995) to estimate a Doppler factor and find T> = 5 (cf. Takahara 1994). Using 
this estimate and their time lag measurements, they derive a magnetic field of B = 0.2 G. From 
non-simultaneous data (Shrader <fc Wehrle 1997), the synchrotron portion of the SED of Mrk 421 
peaks at about ~ 130 eV. Using 

e P eak = 7^-7 2 ^ (54) 

-Dcrit 

, and 7 = (m p /m e )T, we find T ~ 60 and an observer angle ~ 5°. We also obtain an 
equipartition field strength of B eq ~ 10n^ 2 G implying an equipartition parameter of ~ 10~ 2 . 
Although the above value for the bulk Lorentz factor is substantially larger than the mean value 
of (r) ~ 10 found by applying the beaming model to a sample of radio-loud objects (Urry & 
Padovani 1995), its large value may indicate the special nature of gamma-ray loud blazars which 
are characterized not only by small observing angles but also by larger than typical bulk Lorentz 
factors. Despite the crudeness of this calculation, it illustrates the potential applicability of this 
sort of analysis to blazars as well as bursts. 

We thank Jeff Skibo for many useful discussions, and we thank Markus Bottcher for providing 
us with his computer code for calculating the synchrotron emissivity. This work was performed 
while J.C. held a National Research Council-NRL Research Associateship and was supported by the 
Office of Naval Research and the Compton Gamma-Ray Observatory Guest Investigator program. 
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Fig. 1. — Comparison with the analytic result of Dermer & Chiang 1998. The parameters are 
Tq = 300, B = 1 G, n ex t = 1 cm -3 , s = 2. Upper panel: The spectra, from top to bottom, are at 
observer times t Q b s = 1, 10, 10 2 , . . . , 10 6 seconds, with the present calculation plotted as the solid 
curves and the analytic estimate as the dashed curves. Lower panel: Light curves at radio (2.4 GHz, 
dotted), optical (1 eV, solid), X-ray (1 keV, dashed) and gamma-ray (1 MeV, dot-dashed) energies. 
The results from the analytic estimate are over-plotted as symbols. Light curves for each successive 
energy band are displaced vertically by 2 units for clarity. 



Fig. 2. — Upper panel: Radio surface brightness versus perpendicular distance as a function of 
observer angle 6 at constant time delays, £ b s = 1,10, 10 2 ,... seconds (top to bottom). Lower 
panel: Apparent size versus time delay. The data points are from the above calculation and the 
solid curve is equation 39 assuming £ = 3, i.e., a radiative blast-wave. 



Fig. 3. — Numerical results from our blast-wave calculation using model parameters for a burst 
in equipartition, i.e., = 1. Other relevant parameters are £ e = 0.5, Tq = 300, Eq = 10 52 ergs, 
^ext = 1 cm -3 , s = 3, and r\ = 1. Panel a: The bulk Lorentz factor of the blast-wave shell as a 
function of shell radius. The vertical dotted lines indicate the radii of the nearest surface of the 
shell at observer times i b s = 1,10, 10 2 ,... s. The power-law of the bulk Lorentz factor during 
the deceleration phase is C — 1-9. Panel b: Synchrotron and synchrotron self-Compton spectra 
at the above observer times. The spectrum corresponding to each time can be identified by the 
energy, e pea k, of its synchrotron peak which decreases monotonically with time. As in figure [j], the 
vertical lines denote various observer energies — the dot-dot-dot-dashed line corresponds to 1 TeV. 
Panel c: Co-moving electron distributions at the leading surface of the shell at the above observer 
times. Each distribution can be identified by the high energy break which corresponds to the 
synchrotron peak and which also decreases monotonically with time. The dashed curves show our 
analytic estimate (eq. ^2|) for these distributions. Panel d: Light curves for radio, optical, X-ray, 
gamma-ray and TeV energies. 



Fig. 4. — Burst spectra near e pea k for different values of the magnetic field equipartition parameter 
illustrating the effect of different degrees of synchrotron cooling on the burst spectral shape. In 
all three cases, the initial bulk Lorentz factor Tq = 260 (eq. pi]). For the solid, dashed and dotted 
curves we us an electron injection index s = 4. Solid curve: The spectrum obtained for the limiting 
value £b = 5.5 x 10~ 5 (eq. [52] ) which marks the transition between complete and incomplete 
electron cooling. The energy spectral indices (L u oc v~ a ) above and below e pea k are a = —1/3 and 
a = s/2 = 2 respectively. Dashed curve: Here £b = 1 x 10~ 6 , and we have a = (s — l)/2 = 3/2 
above e pea k- Dotted curve: = 1 x 10~ 2 . Above e pea k, a = s/2 = 2, and below e pea k there is an 
additional spectral break due to the cooled electrons which have been injected at the lower cut-off 
energy (eq. |4^) — just below e pea k, a = 1/2 and below the additional break, a = —1/3. Dot-dashed 
curve: £b — 10~ 5 and s = 2. Just above the low energy break we have a = (s — l)/2 = 1/2, and 
the vL v peak corresponds to emission from the electrons at the high energy cut-off. 
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Fig. 5. — Burst calculation with Tq = 260 and £b = 5.5 x 10 5 . All other parameters are the same 
as for figure H The thin curve in panel a is the initial estimate of the bulk Lorentz factor using 
eq. The thick curve is the solution to which our iterative procedure converges. The difference 
in these curves shows the effect of the electrons not being radiative due to the weak magnetic fields. 
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Fig. 3.- 
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Fig. 4.— 
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Fig. 5.- 



